Figure 1(a) presents the top and side views of the geometric structure of the Zn₃Si₂ monolayer, containing three zinc (Zn) atoms and two silicon (Si) atoms in a unit cell. Obviously, the Zn₃Si₂ monolayer has a planar honeycomb lattice structure with P6/mmm (No. 16) symmetry. The black dotted line indicates the minimum repeating unit of Zn₃Si₂ with a = b = 8.05 Å. As shown in Fig. 1(a), there are one unique Zn atom (site symmetry 3g) and one independent Si atom (site symmetry 2d) in the primitive unit of the 2D Zn₃Si₂ monolayer which is structurally similar to the graphene-like honeycomb Be₃Si₂^[39] monolayer. In this planar monolayer, the bond angle Si–Zn–Si is 120° and the Zn–Si bond length is 2.31 Å (as shown in Fig. 1(b)). All Si atoms are trigonally coordinated to surrounding the Zn atoms through sp²-hybridization. Such structure can be seen as transition metal atoms been linearly interposed between two Si atoms in an enlarged silicon framework, it is a benefit to the expansion of the π bond and the stability of the 2D sheet.

Fig. 1.

	Figure Option ViewDownloadNew Window
	Fig. 1. (a) Optimized geometry of Zn3Si2 monolayer, with a unit cell labeled by the black dotted line. (b) A zoom-in figure of (a).

To examined the stability of the Zn₃Si₂ lattice, the cohesive energy^[40] with respect to the isolated atoms is first considered, which is obtained as E_coh = (3E_Zn + 2E_Si – E_Zn3Si2)/5 = 1.95 eV per atom, where E_Zn3Si2, E_Zn, and E_Si are the total energy of the monolayer, the energy of a single zinc atom, and the energy of a silicon atom, respectively. The cohesive energy in our result indicates that the Zn₃Si₂ monolayer might be synthesized by epitaxial growth on an appropriate substrate such as silicone, because the value is comparable with that of grey antimonene (2.30 eV/atom)^[41] and metallic bismuthene (2.26 eV/atom).^[42] In addition, the formation enthalpy also has been considered to investigate the stability of the Zn₃Si₂ lattice further according to the formula ΔH(Zn₃Si₂) = E(Zn₃Si₂) – [ 1 – x]E(Zn) + xE(Si),^[43,44] where E(Zn₃Si₂) is the total energy of Zn₃Si₂ phase per atom and obtained by the first-principles calculations. E(Zn) and E(Si) are the energies of the pure elements Zn and Si, respectively. x is the composition of Si in the compound, . It is found ΔH = –7.3 eV, which implies that the thermal stability of Zn₃Si₂ is good.

The dynamical properties of the 2D Zn₃Si₂ monolayer have been studied by phonon dispersion calculations, as shown in Fig. 2(a). It clearly shows that the 2D Zn₃Si₂ monolayer has no imaginary frequency, indicating that the 2D Zn₃Si₂ monolayer is kinetically stable. Additionally, we also perform AIMD simulations using a supercell of 3 × 3 unit cells (see Fig. 2(b)). Note that the structure still maintains the 2D lattice shape throughout a 3 ps MD simulation at 300 K, indicating that the Zn₃Si₂ monolayer is dynamically and thermally stable at room temperature. Considering that many 2D materials such as graphene, MoS₂, BN, and stanene have been experimentally synthesized, it is expected that a similar method can be used to synthesize a Zn₃Si₂ monolayer.

Fig. 2.

	Figure Option ViewDownloadNew Window
	Fig. 2. (a) Phonon dispersion of Zn3Si2 monolayer, where no soft mode is found. (b) The total energy for the Zn3Si2 lattice as a function of simulation time at 300 K. The inset illustrates the snapshots of the optimized crystal structures of the Zn3Si2 lattice at 1 ps, 2 ps, and 3 ps.

The chemical bonding in the monolayer can be explained by the electron localization function (ELF),^[45] the charge density difference, and the partial charge densities at the Fermi level for this material (as shown in Fig. 3). The charge density difference in Fig. 3(a) shows the hole density near the zinc positions, while the electron is abundant around the in-plane orbitals of silicone, especially along the Zn–Si bonds. The electron density responsible for the formation of bonds is shifted toward the Si atoms. The shape of the distribution in the Zn₃Si₂ monolayer implies that the in-plane orbitals participate in σ-bond formation. The ELF diagram (Fig. 3(b)) shows that the electrons in Zn₃Si₂ monolayer are localized at the Si atoms, whereas the Zn atoms reveal electrons deficiency. The form of electron localization shows the sp² hybridization of the silicone atomic orbitals, which participates in the construction of the σ bonds. Moreover, a small part of the electron density is located in the gap region – TM donates electrons to the metalloid and then a portion of the charge is transferred from the metalloid to the gap region – revealing the characteristics of metalloid bonding. In real space, the partial charge density of the valence band (VB) maximum (Fig. 3(c)) and conductive band (CB) minimum (Fig. 3(d)) confirms that the Zn p_z, d_{xz, yz}, and Si p_z orbitals participate in the formation of the Dirac point. It can further stabilize the 2D framework structure of Zn₃Si₂.

Fig. 3.

	Figure Option ViewDownloadNew Window
	Fig. 3. (a) The charge density difference with an isovalue of 0.02 e/Å−3. The blue and orange mean electron accumulation region and electron depletion region, respectively. (b) ELF of Zn3Si2 monolayer with an isovalue of 0.03 e/Å−3. (c), (d) Top and side views of VB maximum and CB minimum charge density contours at the Dirac point of Zn3Si2 monolayer, respectively.

To evaluate the electronic properties of the predicted monolayer Zn₃Si₂, the projected band structure and the corresponding partial density of states (PDOS) are investigated in the absence of SOC, as shown in Figs. 4(a) and 4(b). The Zn₃Si₂ monolayer is a zero-band-gap semi-metal with the VB and CB touching each other at the K point (Fig. 4(a)). The yellow, gray, blue, and green dots represent the total p_z, p_y, p_x, and d orbitals, respectively. It can be clearly seen that the Dirac point near the Fermi level is mainly contributed by the p_z orbitals. The other two crossing bands are mainly composed of the p_y orbital, which are about 0.5–3.5 eV below the Fermi level. Considering the D_6h point group symmetry of the Zn₃Si₂ monolayer, the Zn-d orbital can be divided into three categories: d_{xy,x² – y²}, d_{xz, yz}, and d_z². In Fig. 4(b) we plot the PDOS of the Si and Zn atoms, respectively. Near the Fermi levels, the VB and CB are mainly from the Si p_z orbital and the Zn p_x,_y, d_{xz, yz} orbitals, respectively. Meanwhile, the VB-1 is primarily derived from the Si p_x,y and Zn p_x,y, d_{xy, x² – y²} orbitals. So, it is considered that the Dirac cone is derived from the special hexagonal lattice structure. As the Zn₃Si₂ monolayer has a planar structure, it is considered that p_x,y and d_{xy, x² – y²} are hybridized to form an σ bond, while p_z and d_{xz, yz} are hybridized to form a π bond. The three-dimensional (3D) valence and conduction bands are also presented in Fig. 4(c), which clearly shows the features of the Dirac cone at its high symmetry K point. The Dirac cone-like electronic states generally mean excellent electronic transport properties. We have calculated the Fermi velocity (V_F) of the Zn₃Si₂ monolayer near the Dirac point by using the equation V_F = ∂E/(ℏ ∂ k). The calculated V_F values along the directions Γ → K and K → M are 2.85 × 10⁵ m⋅s⁻¹ and 1.65 × 10⁵ m⋅ s⁻¹ (without SOC calculation), respectively, which are in the same order of those of graphene (9.5 × 10⁵ m⋅s⁻¹ and 8.2 × 10⁵ m⋅s⁻¹). The corresponding V_F values with SOC are 2.79 × 10⁵ m⋅s⁻¹ and 1.49 × 10⁵ m⋅s⁻¹, respectively, which are similar to the results without SOC.

Fig. 4.

	Figure Option ViewDownloadNew Window
	Fig. 4. The electronic properties of Zn3Si2 without SOC. (a) The orbital-resolved band structures for Zn3Si2 monolayer (yellow: total pz orbitals; light grey: total py orbitals; blue: total px orbitals; olive: total d orbitals). The Fermi level is set to zero. (b) Projected density of states of Zn3Si2 monolayer. (c) 3D projection of Dirac cone in the vicinity of the K point.

The SOC results in “frivolous” splitting of orbital energy levels caused by the interaction of particles’ spin and orbital momentum. The band structure calculated with SOC is shown in Fig. 5(a) It can be seen that the Dirac cone remains. In Fig. 5(b), the zoom-in around Fermi energy shows that a tiny band gap of 1.2 meV appears after considering SOC. Although the crystal lattice of silicene is smaller than that of Zn₃Si₂, the bandgap of silicene (1.55 meV) with SOC is slightly larger than that of Zn₃Si₂ due to 2D silicene materials being undulating and the in-plane p_x and p_y orbitals enhancing the SOC effect.

Fig. 5.

	Figure Option ViewDownloadNew Window
	Fig. 5. (a) The calculated band structures of Zn3Si2 monolayer with SOC. The Fermi level is set to zero. (b) The zoom-in around the Fermi level corresponding to the red dotted box in (a), where the inset indicates the first Brillouin zone.

According to the previous studies,^[46] carrier doping is used as an effective means to manipulate the electronic structure and magnetic properties. The effect of carrier doping on the Zn₃Si₂ monolayer has been studied. It can be seen from Fig. 6(a) that the highest carrier doping concentration is 0.3 electron/hole per atom (carrier density 1.5 × 10¹⁵ cm⁻²). This concentration is experimentally feasible through the electrolysis gates,^[47] which can modulate the carrier density to 10¹⁵ cm⁻² order. The negative and positive values on the horizontal ordinate are for the hole and electron dopings, respectively. It is noted that there is a transition from spin non-polarization to spin polarization through carrier doping in the Zn₃Si₂ monolayer. It is interesting to note that the Dirac half-metallic property is shown in the case of 0.2 hole/atom doping. The band structure of the system with the doping of 0.2 hole/atom is shown in Fig. 6(b). Two fascinating properties, which are independent of the exchange-correlation functional, are observed: (1) the spin-up channel hosts a Dirac cone state at a high-symmetry K point in the vicinity of the Fermi level; (2) the spin-down channel is a semiconductor with a indirect gap of 0.71 eV (shown by a blue arrow in Fig. 6(b)). It is considered that the Fermi energy shifts below the spin-up band, such that the spin-down band is fully occupied, resulting in a truly single-spin Dirac fermion state with 100% spin-polarized currents. Our results indicate that the total magnetic moment of the Zn₃Si₂ monolayer is 1.009 μ_B. However, with the same electron doping concentration, the band splitting of the Zn₃Si₂ single layer is not very obvious (Fig. 6(c)), and the corresponding DHM is not obtained. It means that the spin-polarized Dirac electrons can be obtained by the hole doping even if there is no external magnetic field, which is a unique magnetic phenomenon of Dirac bands.

Fig. 6.

	Figure Option ViewDownloadNew Window
	Fig. 6. (a) The ΔE = E(spin-polarized) – E(spin-nonpolarized) with respect to carrier doping calculated in the Zn3Si2 monolayer. The negative and positive values on the horizontal ordinate are for the hole and electron dopings, respectively. (b), (c) Spin-polarized band structures in the absence of SOC with the doping concentration of 0.2 hole per atom and 0.2 electron per atom, respectively.